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Abstract 

A formulation of Skyrme model as an embedded gauge theory with the constraint deformed 
away from the spherical geometry is proposed. The gauge invariant formulation is obtained 
firstly generalizing the intrinsic geometry of the model and then converting the constraint to 
first-class through an iterative Wess-Zumino procedure. The gauge invariant model is quantized 
via Dirac method for first-class system. A perturbative calculation provides new free parameters 
related to deformation that improve the energy spectrum obtained in earlier approaches. 



1 Introduction 



The nonlinear sigma model (NLSM) with a Skyrme termQ has received a great deal of attention 
after the proposal by Adkins et al[0] that this is a theory for weakly interacting mesons in the chiral 
limit resulting from the more fundamental theory for strong interactions, QCD, in the limit when 
the number of colors iV c is taken very large. Using the collective coordinate method for the isospin 
rotation, the authors of |2| have performed the semi-classical quantization of the model with the 
spectrum of static properties being accurate to within 30% of the experimental data, an extraordinary 
result for such simple model. This model, however, presents some well known shortcomings. There is 
no compelling physical reason to include only the Skyrme term and not higher derivatives potentials. 
Also the quantum picture is known to be troubled by operator ordering ambiguities Q and the results 
following different approaches are not in complete agreement. 

Technically the procedure in || has reduced the field theory problem to that of a free point particle 
on a S 3 sphere. This theory is an example of a nonlinearly constrained second-class system||. The 
quantization of such systems has been intensively studied [§, |5j, §, |7|] in various contexts and the 
relevance of this problem for the quantization on curvilinear surfaces is well appreciate, both in the 
path integral and in the canonical approach. There has been different proposals in the sense of 
improving the physical description of the model|§]. A new possible route, advocated in the present 
paper, would be to deform away from the intrinsic spherical symmetry. However there is no conclusive 
studies of nonlinear constraints of arbitrary geometries available in the literature. This gap makes 
it harder to implement this program since the exact shape of the constraint surface able to improve 
accuracy must come out from a numerical adjustment and is expected to be level dependent. All 
this points up to the necessity of a detailed study of quantization over arbitrary constraint surfaces. 

It is the purpose of this work to illuminate the quantization of Skyrme's collective rotational mode 
resulting from quadratic surfaces deviated from the spherical geometry. While most investigations 
are done towards understanding the quantum nature directly from the 2nd-class formulation, we 
reformulate the model as a gauge theory. We provide the gauge invariant reformulation of Skyrme's 
model for a general nonlinearly constrained second-class surface. This is done through the iterative 
constraint conversion scheme developed by us a few years ago§|. 

The technique that converses second-class constraints to first-class extending the original phase 



space with some auxiliary variables |T(J, was suggested by Faddeev and Shatashivilli [|ll|] with the 
addition of Wess-Zumino terms (WZ) on the original Hamiltonian. As far as we know the treatment 
of nonlinear systems as gauge theories appears already in Balachandran et al in [12] and in a number 
of papers by Yamawaki et al.[I3|. Our presentation here is inspired by the idea of Kovner and 
Rosenstein (KR)|14j], using an analogy with QED to disclose a symmetry hidden in the NLSM. In || 
we interpreted the KR symmetry as the gauge symmetry of a Wess-Zumino extended theory, with 
the geometrical second-class constraints converted to first-class. The hidden symmetry was shown to 
be a residual symmetry of the WZ orbits. The formalism developed in || was used recently in [15 



to study nonlinear models restricted by a more general class of constraints. It provides a clear-cut 
geometrical interpretation for the Wess-Zumino gauge symmetry in the extended space as well as in 
the original configuration space. 

In Section 2 we consider, with detailed care, the four dimensional SU(2) nonlinear sigma model 
with an stabilizing Skyrme term and review the semi-classical expansion of the collective rota- 
tional mode. Reduction to a nonlinear quantum mechanical model depending explicitly on the 
time-dependent collective variables satisfying a spherical constraint is performed. This will set the 
stage to allow the Skyrme's model constraint to be deformed and then turned into a gauge theory in 
the extended WZ phase space. A systematical treatment of the constraint conversion is developed in 
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the remaining of the paper and is distributed as follows. In section 3 the noninvariant aspect of the 
theory is reviewed and our notation is introduced. We consider a general setting and study the mo- 
tion of a point particle in a N-dimensional Euclidean space moving freely on a nonlinear surface Q(q) 
embedded in the 3? . Special emphasis is given to the symmetric generating matrix that embraces 
the geometric information. In section 4 the second-class non-spherical model discussed in section 3 
is transformed to a simple gauge invariant theory, written as the sum of the original plus the WZ 
term. To allow for a concrete calculation we adopt a perturbative point of view and consider small 
deformations away from the spherical symmetry. This introduces new parameters into the theory 
allowing for a better fit with experimental data. The last Section is reserved to discuss the physical 
meaning of our findings together with our final comments and conclusions. 



2 The Skyrme Model Revisited 

A few decades ago Skyrme proposed to describe baryons as topological solutions of the NLSM with 
an appropriate stabilizing term. The semi-classical quantization of the model was obtaining in 
separating the collective coordinate. Let us consider the SU(2) Skyrmion Lagrangian 



L' = I d 3 x \ ^-Tr(dM^U) + ^Tr lu^d^U, U^d u U 
4 v ' 32e z 1 
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where /„- is the pion decay constant and e is a dimensionless parameter. U is a SU(2) matrix 
transforming as U — > AU B~ l under chiral SU(2) x SU (2), satisfying the boundary condition lim U = 
I so that the pion field vanishes as r goes to infinity. There are soliton solutions described by the 
action ([!]) whose topological number are identified with the baryon number. To describe the static 
soliton we start with the ansatz U(r) = exp{if a .x a f(r)} where f a are Pauli matrices, x = x/r 
and lim f(r) = and f(0) = 7T. Performing the collective semi-classical expansion in (P, where 

U(r,t) = A(t)U(r)A^(t) and A e 577(2), we obtain after performing the space integral @|, 



L' = -M + lTr (d Q Ad A- 1 ) . 



(2) 



M and X are the soliton mass and the moment of inertia respectively which, in the hedgehog ansatz 
are given by 
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(3) 



and 




1 = — / dr r sin / 

3 Jo ' " ( - \ \ dr I /'- 

The matrix A may be represented by A = + i a.f satisfying the spherical constraint 

3 



(I; 



1. 



8=0 



In terms of these variables the Skyrmion Lagrangian (0) becomes 



(4) 



(5) 
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3 

l! = -M + 2 1 h l + Afi > ( 6 ) 

with the spherical constraint (d) being implemented by the Lagrange multiplier A. The Hamiltonian 
corresponding to (|6]) is 

H' = M + ^j^tf - \Q, (7) 

82 i=0 

with the canonical momenta, 

IT, = 4Xd,. (8) 
For simplicity we rescale the coordinates and introduce a new Lagrangian and Hamiltonian as 

L = M + L' = \q 2 i+ Afi, 

H = - M + H' = ^-\n, (9) 

where 



Pi 




<b 




A 


= 4JA, 


n 


= Q, 2 - c; 



c = 4J. (10) 

We have reduced the Skyrmion problem to that of a non-relativistic unity mass particle con- 
strained over a S 3 sphere, a well known second-class problem. We are now ready to address the 
question of constraint conversion. A general construction based on the iterative formalism is pre- 
sented in this paper. A particular emphasis is given to the extension of the spherical symmetry with 
dramatic consequences for the energy spectrum. 

3 The Non-Spherical Second-Class Aspect 

In this section we review the main features of a theory for a particle constrained to move on a 2nd- 
class nonlinear constraint surface. Let us consider the mechanical O(N) model which is a theory for 
a point particle with coordinates = (qi, q 2 , ^3, ■ • • , (/tv) moving freely over a nonlinear surface, 

Q = ^qkTk m q m - c = 0. (11) 

Here the surface generator symmetric matrix T^ m and the yet arbitrary constant c characterize the 
surface. We are restricting our study to surfaces bilinear in the coordinates, whose exact structure is 
encoded in T^ m . As so these matrices have a decisive role in determining the constraint structure of 
the model. This role will become clearer as we progress. The classical trajectory and its dynamics 
is governed by the Lagrangian, 
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C= l -q 2 -\^-qTq- C y (12) 

with A being a Lagrange multiplier enforcing (O) as a constraint. To simplify the notation we omit 
the coordinate indices from now on unless to avoid confusion. The canonical analysis gives, 

Pk = Qk, (13) 

and the primary constraint 

0i = Px « 0. (14) 
The corresponding canonical Hamiltonian reads, 

H= l -pl + \^-qTq-cy (15) 

and the primary Hamiltonian is obtained from ([To] ) enforcing the constraint ([14]) with a multiplier 

w, 

H P = H + U(j) 1 . (16) 
The nonlinear surface Q(q) reappears in the Dirac consistency chain as a secondary constraint, 

02 = {01, Hp} 

= a (17) 
The consistency algorithm demands the presence of a tertiary constraint, 

03 = pTq (18) 

and a quaternary constraint as, 

04 = pTp — XqT 2 q. (19) 

where T| m = T kn T nm . Constraints 2 and 3 have a clear geometrical meaning. The secondary 
constraint enforces the particle coordinates q k to take values over the nonlinear surface. The tertiary 
constraint, called as transverse, enforces the particle velocity to take values over the tangent space of 
Q(q) at q k . To disclose the meaning of 0i and 04 we compute the Poisson algebra that follows from 
the complete constraint set, 

{01, 04 } = qT 2 q, 
{02,03} = qT 2 q } 

{0 2 ,0 4 } = 2qT 2 p, (20) 
{03,04} = 2pT 2 p + 2\qT 3 q. 

There are no more constraints. From the time evolution consistency for the 04 constraint, 

0^(qT 2 q)u + {(j )4 ,n} 1 (21) 
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the multiplier u is determined and the consistency chain stops. We may use 04 ~ to determine the 
value of the multiplier A and then compute the Dirac brackets, 



where, 



9km 
hkm 



Here the vector, 



and 



Uk , Qm}* = 
{<?£: , Pm}* = 9km 
{Pk , Pm}* = hm 



(22) 



5km ~ hA-H m 

t^A d m t n p n t m A dkt n p n . 



tk — dk^i 



(23) 



(24) 



A — tktk 



(25) 



are the basic geometric elements of the theory. The constraints 0i and 04 have no physical con- 
sequence. Their presence becomes necessary just to eliminate the multiplier sector and to give 
consistence to Dirac's algorithm. In fact, an alternative approach where this sector is eliminate from 
the outset may be used[0, O. This completes the analysis of the 2nd-class aspect of the model. 



4 Quantization of the Non-spherical first-class Skyrme Model 



In this section we present our main result, i.e., the quantization of the deformed gauge invariant 
Skyrme model and the computation of the energy spectrum. The conversion of the second-class 



constraints to first-class by extension of the phase space |10[ was originally introduced to avoid the 
difficulties involving anomalies in chiral gauge theories! 11], by removing the dynamical degree of 



freedom obstructing gauge symmetry. This scheme has also been used to covariantize the chiral boson 
leading to a system with an infinite chain of lst-class constraints. This formulation 

and dualities 1001 . The 



constraint 18 



has been of some use in recent developments in the study of superstrings| 
logical reasoning behind this approach is to maintain the iterative structure already present in Dirac's 
formalism. In this sense, the iterative method [|T~0| , 15 1 converts the second-class constraints at their 
prompt appearance in the consistency chain. 

To quantize the deformed Skyrme model and then compute the energy spectrum, it is necessary to 
know the non-spherical eigenf unctions corresponding to the Hamiltonian operator that now satisfies 
the new arbitrary geometry. Since it might become so complex, we consider to analyse the Skyrme 
model perturbatively deviated from the spherical geometry, using the spherical eigenf unctions as 
a zeroth-order input without any additional incoveniences. To this end we start proposing the 
perturbative non-spherical Skyrme model as, 



1 11 



(26) 
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where e is a small parameter and the surface generator symmetric matrix Ty is assumed to be 

Tij = Sij + eAij, (27) 

with matrix A^- embracing the deviation from the spherical surface. The Hamiltonian computed 
through Legendre transformation is, 

U = \pl + A (^i - -y 'l.Xj'U ~ c). (28) 
As in section 2, the model has four second-class constraints, 



Pa, 

1 2 1 A 

q.p + eq.A.p, 



(29) 



P 



Recall that the first and the last constraints have no geometrical meaning. They are in fact artificial 
relations generated by the Dirac iterative process to preserve the original consistency presente in the 
dynamical phase space. Hence, these four second-class constraints could be reduced to only two, the 
geometrical ones (02,03)- Computing the Lagrange multiplier A from the last constraint in (|29|), to 
first-order in e gives, 



Bringing back this result into the Hamiltonian fl2q), it becomes 



(30) 



n 



2 2 , 

QjPi + £ 



4c 



p_ 

4c 



(g.A.g) 



p p 1 s 

— (g.A.g)g| + —(g.A.g) - -p.A.p + —(p.A.p)q 



(31) 



with two second-class constraints, denoted subsequently as Xi — 02 and \2 = 03- At this stage the 
gauge invariant formulation of the second-class model with the introduction of the canonical WZ 
variables (9, Kg) starts|10|, Firstly changing the nature of the geometrical constraint from second 
to first-class, 



Xi -> Xi 

X2 ~> X2 



Xi + 9 

X2 - Ke, 



(32) 



so that the Poisson bracket between them becomes, 



(33) 



After that, it is imperative to compute the new Hamiltonian that preserves the Dirac iterative 
process in the extended phase space, 



{xi,H} 



a E 



X2- 



(34) 
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Therefrom the Hamiltonian counter-term is found and, subsequently, the gauge invariant Hamilto- 
nian, 



H=?f + e 
Ac 



P 



P 



Ac 



j-(q.A.q) - ^(q.A.q)q? + — (q.A.q) - -p.A.p+—(p.A.p)q* 



Ac 



--r(q 2 + e{q.A.q)){q 2 + 2e{q.A.q)) 
Ac 



(35) 



This concludes the gauge reformulation of the deformed Skyrme model. 

To conclude our task, this gauge invariant model has to be adequately quantized. Thereby, the 
WZ symmetry is partially fixed as, 



i> = e. 

Consequently, x 2 turns to second-class constraint, 

{0,X2} = q 2 + 2eq.A.q 
that allows to compute the corresponding Dirac brackets, 



(36) 



(37) 



{qiAjY = {Pi,Pj}* = Q, 

UhPj}* = hj- (3* 



Notice that these brackets turn out to be canonical, illustrating the Maskawa-Nakajima theorem pT 
The reduced Hamiltonian is rewritten as, 

^ ~ M _ , _ I A N P 2 I A \ 2 , Q 2 , A \ 1 A / ^1- A -P , (q-PY 



H = —piMijPj + e 



(q.A.q) - —(q.A.q)q + —(p.A.p) - -p.A.p - (q.p) — + — - (q.A.q) 



2c ' 8c 2 v ' Ac ' T 2c 2cq 

(39) 

where the matrix M^-, 

M 4 , = S V - "-f (40) 

is singular. Xi|e=o = Xi remains a first-class constraint, identified as the "Gauss Law" generator of 
the remaining symmetry. 

At this point we are interested to compute the energy spectrum of the model from the invariant 
Hamiltonian above and to compare it with the one given by ANW@. To calculate the energy 
spectrum we will use the Dirac quantization method for first-class constraint systems. To this end 
we impose the first-class constraint in the spatial sector of the Hilbert space, 

(q 2 - c)\phys >= 0. (41) 

The momenta representation must reflect the constraint presence and its associated algebra. The 
first-class nature of the geometric constraint is represented as, 

[q 2 ,p i ] = 2ihq i . (42) 
This algebra is satisfied by the following commutator, 



[QhPj] = iftdij, 

with the canonical representation for the momenta operator p i: 



(43) 



Pi = -ihdi, (44) 

which satisfies the relations ( ^ - |43| ). The spectrum of the theory is obtained by evaluating the 
Hamiltonian meam value with typical eigenfunctions defined, for instance, \polyn >= -^-W (q'i + ^2 ) z @] 
to give, 

H\polyn >= Ei\polyn > . (45) 
A direct calculation produces in this case 



H\polyn >= —1(1 + 2) 



q.A.q 



h + iqz) 1 + 2e(gi + iq 2 ) 1 1 q b A bc d c (q 1 + iq 2 ) 



Therefore, the consistence between ( [4"5]) and fl46| ) requires, 

(Qi + iq2) l ~ 1 qbA bc d c (qi + iq 2 ) = {qi + iq 2 ) 1 ■ 
For instance, a simple choice for the matrix A cb , turning the surface slightly oblate, 



(46) 



(47) 
















1 














1 





\0 








0/ 



(48) 



produces the following eigenvalue, 



E, 



4^ + 2) 



1 + 2e 1 



2c 



(49) 



Clearly, the extra term reflects the correction generated by the deformation away from the spherical 
symmetry. Note that this deformation matrix Aij has produced a new free parameter that may be 
used to better the fit energy spectrum to experimental data. Notice that for e > the above result 
improves the ANW estimate 



Ei = E. 



ANW 
I 



+ eAEi 



(50) 



Before closing this section we would like to mention an important feature related to the inter- 
nal consistency of the formalism. If we choose the deformation matrix function Ajj = 5{j (which 



automatically satisfy ([47])) then such deformation corresponds only to a change in the radius of the 
sphere. A direct calculation from ( [46]) shows that the spectrum remains unchanged. This is im- 
portant since the radial direction was seen to be the gauge orbit [T^j - spheres of different radii are 



physically equivalent. Only deformations on the geometry produce energy corrections. 
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5 Final Discussions 



The quantization of the Skyrme model is a well known example of quantum mechanics on curved 
space. This comes from the presence of nonlinear second-class constraints but the passage to the 
quantum world is plagued by ambiguities. In general the use of the Poisson bracket structure is in 
conflict with the constraints. The use of the Dirac brackets helps solve this problem classically, but 
the new algebraic structure ends up being coordinate dependent posing a new problem in terms of 
ordering ambiguities at the quantum level. Different representations for the momentum operator 
brings additional terms into the quantum Hamiltonian associated to the zero-point motion. The 
quantum dynamics becomes dependent on the particular ordering adopted and is not unique. 

To cure these sort of problems we propose to treat nonlinear second-class systems as gauge theories 
by converting the constraints into lst-class. After conversion has been accomplished, the quantization 
may be implemented by the Dirac or Gupta-Bleuler procedures where the lst-class gauge constraints 
select the physical Hilbert space. Since the Poisson brackets remain as the underlying algebraic 
structure, no ordering ambiguities affect the quantization process. The question that remains is 
how to convert efficiently the constraints so that different geometries could also be treated. This 
seems important to improve the quantization of the rotational modes of the Skyrme theory. The 
development of a formalism to answer this question is the main contribution of this work. 

The development of this formalism was done in Section 4. We studied the model with the most 
general quadratic constraint and showed that a single WZ term produces the desired covariantization 
of the Skyrme Lagrangian, even for arbitrary geometries. The exact computation of the spectrum 
for an arbitrary geometry has shown to be a too difficult problem. We proposed to handle it in 
a systematic fashion by perturbative technique. A variety of deformations may then be examined 
through a proper choice of the matrix Ay satisfying the condition fliTD - The new parameters related 
to the geometry deformation modify the energy spectrum. Notice that when these parameters pre- 
serve the spherical symmetry the old results of ANW are recovered. It should be noticed that this 
deformation occurs in the internal space of constraints and, therefore, can not be compared directly 
with the results of the "deformed skyrmions" proposed earlier where the deformation was encoded in 
the profile function of the hedgehodge solution || Our solution provides new free parameters 

related to the deformation which may be adjusted adequately to fit experimental results. 
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